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Abstract 

The first and second Born approximation are studied with the path integral rep- 
resentation for T matrix. The T matrix is calculated for Woods-Saxon potential 
scattering. To make corresponding integrals solvable analytically, an approximate 
function for the Woods-Saxon potential is used. Finally it shown that the Born se- 
ries is converge at high energies and orders higher than two in Born approximation 
series can be neglected. 

Keywords:path integral; scattering; S matrix; Born approximation; Woods-Saxon 
potential 

PACS numbers: 31.15.xk; 13.85.Dz; 11.55.Jy 



* Email: m.pahlavani@umz.ac.ir 
^Email: r.morad@umz.ac.ir 



1 



1 INTRODUCTION 



There is a general idea that there are several ways to describe nature. In spite of equality 
of these ways, they are different in prediction of new laws of physics. Before testing 
the power of the predictability of a new formalism, in the first step, it is necessary to 
test it with other laws of physics. We are all familiar with the standard formulations of 
quantum mechanics, developed more or less concurrently by Schroedinger, Heisenberg and 
others in the 1920's, and shown to be equivalent to one another soon thereafter. In 1933, 
Dirac made the observation that the action plays a central role in classical mechanics (he 
considered the Lagrangian formulation of classical mechanics to be more fundamental than 
the Hamiltonian one), but that it seemed to have no important role in quantum mechanics 
as it was known at the time. He speculated on how this situation might be rectified, and 
he arrived at the conclusion that (in more modern language) the propagator in quantum 
mechanics corresponds to exp(iS/h), where S is the classical action evaluated along the 
classical path. In 1948, Feynman developed Dirac's suggestion, and succeeded in deriving 
a third formulation of quantum mechanics, based on the fact that the propagator can be 
written as a sum over all possible paths (not just the classical one) between the initial 
and final points. Each path contributes exp(iS/h) to the propagator. So while Dirac 
considered only the classical path, Feynman showed that all paths contribute: in a sense, 
the quantum particle takes all paths, and the amplitudes for each path add according 
to the usual quantum mechanical rule for combining amplitudesfl, 2]. Path integrals 
give us no dramatic new results in the quantum mechanics of a single particle. Indeed, 
most if not all calculations in quantum mechanics which can be done by path integrals 
can be done with considerably greater ease using the standard formulations of quantum 
mechanics, path integrals turn out to be considerably more useful in more complicated 
situations, such as field theory. But even if this were not the case, it is believed today 
that path integrals would be a very worthwhile contribution to our understanding of 
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quantum mechanics. Firstly, they provide a physically extremely appealing and intuitive 
way of viewing quantum mechanics: anyone who can understand Youngs double slit 
experiment in optics should be able to understand the underlying ideas behind path 
integrals. Secondly, the classical limit of quantum mechanics can be understood in a 
particularly clean way via path integrals. It is in quantum field theory, both relativistic 
and nonrelativistic, that path integrals (functional integrals is a more accurate term) 
play a much more important role, for several reasons. They provide a relatively easy 
road to quantization and to expressions for Greens functions[3, 4, 5], which are closely 
related to amplitudes for physical processes such as scattering and decays of particles. 
Furthermore, the close relation between statistical mechanics and quantum mechanics, 
or statistical field theory and quantum field theory, is plainly visible via path integrals. 
Now, path-integral formalism is widely used in many branches of theoretical physics and in 
particular in nuclear physics[6]. The existence of an strong nuclear force in atomic nuclei 
revealed by the exceptional role of the nuclear magic number provides the foundation of 
the nuclear shell model. This strong force is believed to be approximated most closely by 
a Woods-Saxon potential [7] either from analyzing the radial dependence of the nuclear 
central force or by deriving it from a microscopic two body force acted in neutron proton 
scattering [8]. The spherical Woods-Saxon potential that was used as a major part of 
nuclear shell model, was successful to deduce the nuclear energy levels [9]. Also it was 
used as central part for the interaction of neutron with heavy nucleus [10]. The Woods- 
Saxon potential was used as a part of optical model in elastic scattering of some ions with 
heavy target in low range of energies [11]. In this paper we will use the path integral 
representation of the T matrix in potential scattering for driving the different orders 
of Born approximation with the Woods-Saxon potential. This representation is not a 
phase-space path integral but it is a particular path integral over velocities that reduced 
the complexity of path-integral [12]. In the limit of large scattering times where energy 
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conservation, the " dangerous" phases from S matrix are created. The " Phantom" degrees 
of freedom is utilized to eliminate these phases. In addition, energy conservation is applied 
by imposing a Faddeev- Popov-like constraint in the velocity path integral [12]. In Sec. II 
we outline the T matrix representation from the velocity path integrals. In order to 
obtain different terms of Born approximation, we use path integral representation for the 
scattering in presence of Woods-Saxon potential in Sec. III. In this investigation we show 
that the Born approximation is a converge series, therefore in elastic scattering analysis, 
we can neglect the orders higher than two in Born series. 



2 Path Integrals for the T matrix 

In the framework of nonrelativistic potential scattering, consider a central potential, V(r) 
that vanishes at infinity. k« and k/- are the initial and final momentum of a particle with 
mass m (see Figure. 1). In this calculation, supposed the scattering states are normalized 
and h = 1[12], so 

(<p s \tp i ) = {2«)*P{k i -k t ). (1) 

The S matrix is the matrix element of the evolution operator in the interaction picture 
that taken between scattering states and calculated at asymptotic times: 



*/ 



lim 

T — >oo 



Ui(T,-T) 



lim e W%) r 

T — >oo 



U(T, -T) 



(2) 



where the Ui(T, — T) is the time evolution operator that is defined as follows 



Ui(t b , t a ) = e iAotb exp \-iH (t b - t a 



(3) 
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Figurel: Scattering geometry for a potential of radius R, the impact parameter b. Incom- 
ing and outgoing momenta are kjj , and the mean momentum is K = (k» + kf)/2. 



The T matrix is defined usually by subtract the identity from S matrix and factor out 
an energy conserving Dirac delta function 



Si^ f = (2tt) 3 5( 3 ) (k< - k f ) - 2m5 {E { - E f ) %. 



*f ' 



(4) 



where E{ = Ef = E = k 2 /(2m) is the common scattering energy. By integrating func- 
tionally over velocities instead of paths[14,15,16] the following formula is achieved[13] 



i—*f 



lim exp ( i—T\ [ d 3 r e~ iC ^ r J\f 3 (T, —T) ( V 3 v exp i f +T dt—v 

n — >oo y Am J J J J-t 2 

i [ T dtV ( r + — t + Xolt) 
J-t \ m 



x < exp 



where 



J Vv exp \i jf dt-v 2 (t) 



(6) 



and 



Xv(t) = ~ dt'sgn(t - t') w{t') . 



(7) 



The momentum transfer and the mean momentum are defined by 



q = k f - ki , K = - (kj 



(8) 



It can be shown that the 'dangerous phases' q 2 T/4m that are proportional to T, are can- 
celed in each order of perturbation theory by introducing 'Phantom' degrees of freedom( 
dynamical variables with the wrong sight kinetic term) [12], so that the limit T — > oo can 
indeed be applied. Then the following path-integral representation for the S matrix is 
obtained [12], 

(S - ll^ f = lim / d 3 r e - iC L T LA/YT, -T)| 6 / V 3 vV 3 w exp i [ +T dt ^ (v 2 (t) - w 2 (t) 

■I T >oo J J J-T 2 ^ 

•jexp -ij + ^dtV (r + + x v (t) -x w (0)^ -1 

where w(t) is a three-dimensional ' ant i velocity'. This is similar to the Lee- Wick approach 
to quantum Electrodynamics where they introduced the fields with a wrong sign kinetic 
term to remove all infinities[17,18,19]. To extract the T matrix from the S matrix for 
weak interaction, we can develop in powers of the potential and factor out in each order 
an energy conserving 5 function. But to achieve this without a perturbative expansion of 
S matrix usually the trick which Faddeev and Popov have introduced in field theory for 
quantization of non-Abelian gauge theories is used [20] and the following expression for 
the T matrix is obtained[12], 

K 



(3-3) 



m 



d 2 be~ iC L h |A/f / VhiV'w ex V {i I dt 



m r 



v 2 (t)-w 2 (t)} 



(10) 



x ! e ix K (b,v,w) 



Here the limit T — > oo is taken and the corresponding Gaussian normalization factor is 
written as 



M := A/"(+oo, -oo) 



'ID 
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In this equation the phase %k ^ s defined as 

X K (b,v,w) = -/ dtV [b + -t + x v (t)-x w {0)-\K , (12) 

where b is the impact parameter. 9 is the scattering angle,one may obtained 

q= |q| =2fcsin^J , K = |K| = A; cos . (13) 

For t = that is the most symmetric choice the gauge parameter, A = Kto/m, is zero. 
The superscript "3-3" indicates that there are three-dimensional antivelocity that used 
to cancel divergent phases in the limit of asymptotic times moreover the three- dimen- 
sional velocity variables. Then by expanding the exponent in powers of the potential the 
complete Born series is reproduced from this path-Integral representation as follows[12], 



7^/=:£T B , (14) 



n=l 



with 



T » (3 " 3) = £ / d%e ~ H (/_ T * J ffi nr.) ) d5) 



exp(z^p r [b + x^)]) Gi 3 " 3 ). 



8=1 



Where Gl 3 3 ) is evaluated as 



= exp { -J- £ Pi ■ Pj (T - \t t - tj\ - T) \ , (16) 



and 



x rcf (t) = -t. (17) 
m 

V(pj) is the fourier transform of potential 

V(p.) = J d 3 rK(r)e i( l- r . (18) 
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3 Born Approximation 

In this section, we calculated the first and the second order of Born approximation from 
the path integral formula for T matrix. For first Born approximation from the Eq.(15) 
we have 

r/ 3 " 3 ) = K-\ dHe-^ [ +T dt f v-( p ) x exp Up • (b+ -t)\ . (19) 
m J J-t J [2-k) 6 [ \ m ) J 

That with integration over b and t, we obtain 

= V(p ± = q, P|| = 0) = V(q). (20) 

Where p ± and py are the components of p perpendicular and parallel to K respectively. 
Also q is a vector in the plane which is perpendicular to K. The second Born approxi- 
mation is calculated in the following formula 

rt Z) = Z ^ JL I I d 3 Pl d 3 p 2 V( Pl )V(p 2 )5^ (q - p 1± - p 2± ) (21) 
m J J 

x §( K \ Pl ' P2 ) 5 (p — Pl ' p A 

V 1 m 2m J \ 2 m 2m J 

To simplify the above equation, we suppose that p : • p 2 = 0, then the Eq.(21) can be 
rewritten as 

rr 3) = ^P f* P v<p)v(<i-p)- (22) 

Integration over p can been done in the plane perpendicular to K. Now, we apply 
this formalism to calculate the first and second Born approximation for Woods-Saxon 
potential. The Woods-Saxon potential can be defined by 

v <r) = rz V*.v (23) 

1 + exp(^p) 

This potential may be represented, with precision better than 3% for any r value ( see 
figure. 2 ) by 



C(x < 0) = 1 - - e x + - e 2x (25) 
8 8 

C(x>0)=e-*(1 ~le' x + le~ 2 *). (26) 
This approximation is particularly useful in obtaining analytical expression for integrals 




Figure 1: Nuclear potential for the 16 + 20 Ne scattering represented by Woods-Saxon 
potential(solid line) and the approximate function VqC{{t — Ro)/a) (symbol •). 

that involve the Woods-Saxon potential. A such approximation is used for tow parameter 
Fermi distribution density functional]. With this approximation the fourier transform of 
Woods-Saxon potential is 

j p -Ro/a c -2R /a 

"^'^'WTIF-^ 1 (27) 

where q is the magnitude of q. And the second Born approximation can be written as, 

T 2 (3 " 3) = [49e^/^ 1 (g) - A2e~^g 2 (q) + 36e-^g 3 (q)l (28) 

where gi(q), #2(9) and (73(g) are solvable integrals over p, but their results are too compli- 
cated, thus their results are not appeared here. It should be noted that q is a function of 
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scattering energy (E) and angle (0), therefore 



9i(q) = / dp 



P 



(29) 



(p 2 a 2 + l) 2 ((q-p) 2 a 2 + 1) 



2 ' 



92(g) = dpp{ 



1 



1 



], (30) 



(p 2 a 2 + l) 2 ((q - p) 2 a 2 + 4) 



2 + 



(p 2 a 2 + 4) 2 ((g-p) 2 a 2 + l) 2 



93(q) = / ^ 



P 



(31) 



(p 2 a 2 + 4) 2 ((g-p) 2 a 2 + 4) 2 ' 



We wish to calculate the first and second Born approximation for the 16 + 20 Ne. The 
parameters of Woods-Saxon potential for this system are presented in table. 1. The first 
and second Born approximation, for this reaction, as a function of scattering energy and 
angle are shown in Figure. 3 and Figure. 4 respectively. The graphs show that the values 
of Ti and T2 are significant in small angles and extremely reduced with increasing of 
angle. For comparison between first and second Born approximation, we show these 
approximation versus the scattering angle 9, at E=24.5 MeV in Figure. 5. It can be seen 
from this Figure that the value of T2 is smaller than T\ and can be neglected. Also in fig. 6 
the values of T\ and T2 for 9 = 0°, 9 = 45°, 9 = 90° versus scattering energy are plotted. 
These graphs show that the value of T 2 , in all angles and in energies below and near the 
coulomb barrier, is very smaller than T\. This is arises from fact that the Woods-Saxon is 
a short ranged potential. In low energies the repulsive coulomb force prevents the nuclei 
to be close. So it is not convenient to use the nuclear potential to study the scattering. 
In this energies, the scattering is Ratherford scattering. Finally we can conclude that 
the Born approximation is a converge series and we can neglect higher order of Born 
approximation in elastic scattering studies. 
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Table 1: Woods-Saxon parameters for 16 O + 20 Ne reaction. (i? = ^o(^i + ^2 )) 
[22] 



System 




r (/m) a(/m) 


ie +20 ^ 


10 


1.48 0.45 




Figure 2: First Born approximation versus scattering energy and angle(#) for 16 + 20 Ne. 
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Figure 4: The first and second Born approximation calculated with path integral formal- 
ism for 16 + 20 Ne scattering versus 9 in E = 24.5 MeV. solid line shows the first Born 
approximation and dash line shows the second Born approximation. ( in scale 10~ 7 ). 
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Figure 5: The first and second Born approximation calculated with path integral for- 
malism for 16 + 20 Ne scattering versus scattering energy (E) in three angle : 8 = 0°, 
9 = 45° and 9 = 90°. solid line shows the first Born approximation and dash line shows 
the second Born approximation. ( in scale 10~ 7 ). 
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